Mechanism for the a — V e phase transition in iron 
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The mechanism of the a-e transition in iron is reconsidered. A path in the Burgers description of the bcc/hcp 
transition different from those previously considered is proposed. It relies on the assumption that shear and 
shuffle are decoupled and requires some peculiar magnetic order, different from that of a and e phases as 
found in Density-Functional Theory. Finally, we put forward an original mechanism for this transition, based 
on successive shuffle motion of layers, which is akin to a nucleation-propagation process rather than to some 
uniform motion. 
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I. INTRODUCTION 

The bcc-hcp transition in iron has recently been the sub- 
ject of intense experimental 1 ^* and theoretical^— works. At 
room temperature and pressure, a-iron is a bcc metal with fer- 
romagnetic (FM) order. Upon pressure, iron exhibits a phase 
transition at ~ 1 3 GPa to hep structured— with no magnetic 
order. 17 The hysteresis at the transition is large, 18 and as pres- 
sure is changed, the transformation occurs rapidly, suggest- 
ing that the transition is non-diffusive and martensitic. Recent 
work underlines the importance of antiferromagnetic (AFM) 
fluctuations in hep e iron 4 . 

In order to describe a iron from an ab-initio point of view, 
Dynamical Mean Field Theory 1219 in combination with Den- 
sity Functional Theory (DFT) or other advanced schemes 13 
have been used. In particular, the localization of elec- 
trons in the non-bonding e g state but not in the bonding t2 g 
state is a key point 19 , for understanding paramagnetic bcc 
iron. However, DFT ab initio techniques alone can handle 
FM bcc iron and the thermodynamical fundamentals of the 
a — > e transition 20 : The Generalized Gradient Approximation 
(GGA) correctly captures the basic Stoner mechanism for the 
reduction of magnetism, under the increase of atom coordi- 
nation during the a — > e transition, as well as the larger spin 
polarization of e g orbitals with respect to t2 g in the a phase. 

The a — > e transition can be viewed 2 - as combining an 
anisotropic compression (shear) in the (100) direction of bcc 
with a shuffle in the (Oil) plane, that corresponds to a zone- 
boundary transverse phonon mode (see, e.g, Ref. ,5). Due its 
martensitic nature, the transformation is triggered above the 
thermodynamical threshold, i.e., when the a and e enthalpies 
are equal, which appears above ~ lO-OGPa, 5,7 ^ as computed 
in the GGA. 

However, its ab-initio description is debated. Although the 
broad outline of the transition has been settled, a detailed de- 
scription faces the problem of accurately computing the rele- 
vant energy barrier. Ekman et al. 5 demonstrated that the tran- 
sition is first-order, with no dynamical precursors, and that 
it is caused by the effect of pressure on the magnetism of 
iron. Johnson and Carter 7 have shown that a minimization wit. 
shuffle and shear, considered as independent variables, leads 
to a cusp in the energy profile with a particularly low energy 
barrier This pathway imposes a discontinuity of the shuffle 
displacement, at no additional energetic cost. Rejecting this 
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possibility, Liu and Johnson 9 stated that shuffle and shear are 
coupled and cannot be minimized separately. Still, the barrier 
energy found this way is too high — within the range of pres- 
sures for which the transition is expected to occur — to permit 
the transfoimationi& No mechanism allowing for a transfor- 
mation with a reasonable energy barrier has been proposed. 

Investigations have established the role of complex mag- 
netic structures for some transition paths of the a — > 7 (Refs. 
[Tol I22I and I23I) and a — >• e transitions, 6 e.g., for e— Fe spin 
spiral states have the same energy as AFM structures. 8 Yet no 
study has been carried out on the effects of magnetic oideiing 
on the pathway of the bcc-hcp martensitic phase transition. 

As to the kinetics of the transition, no intermediate states 
were detected in experiments on the nanosecond time scale, 24 
which suggests that the transformation is very rapid and prob- 
ably propagates at nearly sound velocity. Moreover, it has 
been shown recently that for cobalt, the transition occurs 
through microstructural avalanches. 25 Such behavior implies 
that the transition is not of the military type but instead in- 
volves localized tiansfoimation events. No such possibility 
has been considered yet for iron. 

In this paper, we first propose (Sec. [TTJ an alternative de- 
scription of the path for the a — > e transition within the Burg- 
ers mechanism. Second, using GGA calculations, we show 
the importance of AFM order to describe the eneigetics of 
this path (Sec. IIII Al ). Thiid, and this is oui main result, we 
emphasize (Sees. iHTBl and fTTTCb that the shuffle mechanism 
could occur layer after layer in a non-simultaneous way, thus, 
bridging the gap between the low-energy pathway of Ref. 7 
and a description of the transition without any discontinuity in 
displacement^ 



II. COMPUTATIONAL SETUP AND METHODOLOGY 

The two-atom unit cell is replicated four times in the 
(OOl)hcp || (HO)bbc direction into an eight-atom supercell, 
with periodic boundary conditions (PBCs). The unit cell is 
represented in Fig. [T]with lattice vectors R\, R2, and R3. 
Three parameters are necessary to determine the unit cell: the 
angle a, the value of c = i?2/4 and the value of R\. Three 
macroscopic quantities can be used to compute them: the vol- 
ume, the c/a — i?2/(4i?3) ratio and the shear e. In our paper, 
the volume, namely, 71.5 bohr 3 /atom, is such that hep and bcc 
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FIG. 1. (Color online) Scheme of the (OOl)hcp || (llO)bbc plane in the 
bcc configuration. The vectors Ri, R2, and .R3 define our simulation 
cell. The red atoms and blue atoms lie in two different planes (see 
Fig. |2j. R2 is perpendicular to .Ri and R3, and R2 = 4c. The shear 



is defined by e = (2 sin a - ^)/(l — -^). 
and e = 0. In hep, 2 since = 1, thus, e — 1 
conservation of volume (R2R1R3 sin a = V) 



In bcc, 2 sin a = -4= , 
From this definition, 



conservation of - 

a 

( = = ^f)' ant ^ me g eometr i c relation 2i?3 cos a = Ri, give 
Ri, R2, and .R3 as a function of the value of the shear. 
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FIG. 2. (Color online) Sketch of possible paths for the bec/hep tran- 
sition in the shear shuffle plane. TS is the transition state, and the 
black line corresponds to the lowest pass to go from the bcc valley to 
the hep valley. The red and green lines are the transition paths that 
are studied in Fig.|4](see also text), ei = 0.5 and £3 = 0.875. 



energies coincide. The c/a ratio is kept constant and equal to 
y/8/3 ~ 1.633, which is the value in the bcc structure and 
in the ideal hep compact structure. The shear is defined by 
e = (2 sin a — -2=)/(l — -2=) and is equal respectively to 
and 1, respectively, in the bcc and the hep structures. For each 
value of the shear, we computed the energy for three values 
of the c/a ratio, namely, 1.605, \/8/3 ~ 1.633, and 1.668 to 
check that the sheared bcc structure is stable with respect to a 
change of c/a and that c/a = y/8/3 is close to the minimum. 
The results are insensitive to the small difference between this 
value and the true minimum. Finally, the shuffle 77 is propor- 
tional to the distance of one blue atom between its position 
during the transition (see Figj2j and the reference position in 
bcc. 

We focus our paper on the shuffle part of the Burgers mech- 
anism. Because of the different physical timescales involved, 7 
we assume that the lattice shear deformation mode (e) and the 
shuffle mode (77) are uncoupled in such manner that, with re- 
spect to shear, shuffle is instantaneous. Volume and c/a ra- 
tio are, thus, kept constant during the shuffle. The path pro- 
posed in Ref. 7 is defined by the equality of the energies of 
the sheared bcc and hep phases. In Ref. 7, however, the en- 
ergetics of the shuffle is not explored. On the other hand, Liu 
and Johnson 9 couple shear and shuffle modes. Their path is 
schematically represented in black in Figj2] In this figure, the 
shear (ei) is fixed to half the value necessary to go from the 
bcc a phase to the hep e phase. It is the shear at the transition 
state computed by Ekman et ali^ at constant volume. This 
choice guarantees that the shuffle mode goes through the tran- 
sition state (TS) of Ref. 5 and corresponds to the energetically 
most favorable path, represented in solid red in Figj2] 

In this paper, we are concerned by the shuffle mechanism 
along the BC line (see letters in Fig. |2). Since the shuffle of 



a large number of atoms can be described in several different 
ways, our aim is to understand whether the shuffle of all the 
layers is simultaneous or not. 



A. Simultaneous shuffle 

Let us first consider simultaneous shuffle (S). By 
AAAAAAAA, we denote the stacking of the eight atoms in 
the bcc sheared phase, and by BBBBBBBB, we denote their 
stacking in the hep sheared phase. In this notation, A and 
A, and B and B, correspond to atoms in the bcc- or hep-like 
configurations, respectively. The atoms in the A and A lay- 
ers have 8 neighbors, whereas atoms in the B and B layer 
have 12 neighbors. During the simultaneous shuffling of all 
four atoms of the A-type layer into a B-type layer, A-type 
layers are transformed into B-type layers. The transformation 
is, thus, AAAAAAAA ->• BBBBBBBB [Fig. Eft)]. 



B. Consecutive shuffle 

Consider next the consecutive individual shuffle of layers 
(C), schematically represented in Fig.|Hb). First, we compute 
the individual shuffling of layer 2 from A to B. It corresponds 
to the transformation AAAAAAAA -> IBIAAAAA. Due to 
the PBCs, this creates intermediate (I) configurations for the 
first and third layers. It corresponds to an atom layer sand- 
wiched between one B and one A layer. Their local environ- 
ment is, thus, neither bcc nor hep. 

In a second step, starting from the latter configuration, we 
study the shuffling of the second layer up to the new config- 
uration IBBBIAAA. During this transformation, the fourth 
layer shuffles from configuration A to configuration B. Thus, 
the third layer now has a full hep-like configuration around it 
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FIG. 3. (Color online) (a) Total energy for the simultaneous shuf- 
fle Si computed for FM, AFMI, AFMIIri and Non-Magnetic (NM) 
configurations. S stands for Simultaneous and the index 1 refers to 
the shear ei (see Fig. [2}. The energy is given for a supercell of eight 
atoms to allow for an easy comparison with Fig. [4] (b) Stacking of 
the layers during the transition. In black, the A and A layers. In 
orange, the B and B layers. The black (respectively, orange) color, 
thus, indicates that the first shell of neighbors is bcc- (respectively, 
hep-) like. In the inset, the evolution of the magnetic moment for the 
eight atoms during the FM — > AFMI transition. 



in agreement with Ref. |g where, however, AFMI order was 
not considered. This confirms the importance of magnetism 
for this transition: 5 ~ 7 ' 9 Without any change in magnetism, the 
transition would not occur. In the rightmost part of the curve, 
the energy decreases to a local minimum corresponding to a 
deformed hep phase. We computed the same curve with a cell 
of four atoms allowing for the creation of AFMII order^ We 
find, in agreement with Refs. |32| and [ft that AFMII is more 
stable for the hep structure (p=4) compared to AFMI. How- 
ever, and interestingly, we also find that past the cusp, the 
AFMI structure is more stable than the AFMII. This highlights 
the important role of complex magnetism for the DFT descrip- 
tion of the transition path and not only for the stable structure. 5 
A complete study of magnetism, especially noncollinear, lies 
beyond the goal of our paper. We nevertheless show that an 
AFMI description of the transition path is particularly adapted 
to the present case. We call E Sl (p) the combination of the 
FM and AFMI curves. The energy profile qualitatively agrees 
with previous results, involving a smaller unit cell with two 
atomic layers (Figs. 5 and 6 of Ref. 5) even if Ekman et al. 
have considered a NM ground state for the hep phase and the 
PW91 GGA functional. 



B. Consecutive shuffle 



and now is, thus, labeled B. The process can be pursued until 
all layers have shuffled. 



C. Computational details 

We use the Projector Wave Augmented 26 implementation^! 
of ABINIT (Ref.lH) in the GGA PBE approximation of DFT. 
Atomic data include 3s and 3p semi-core states,— the cut- 
off radius is 2.0 a.u., the energy cutoff for the plane-wave 
expansion is 20 Ha, and the convergence criterion for the 
charge-density residual is 10 -9 . With these atomic data, re- 
sults of Ref. (271 ar e reproduced. A 18x4x26 special /c -point 
mesh 29 and a Gaussian smearing of 5.10 -4 Ha for electronic 
occupations 30 are used in order to obtain a good estimate of 
the energetic and magnetic properties. 



III. RESULTS 

A. Simultaneous shuffle 

Energies for the simultaneous mechanism are displayed in 
Fig EI a). The reaction coordinate (p) is taken proportional to 
the shuffle (77) of the supercell atoms and varies from (bec- 
like) to 4 (all atoms in hep-like environment). Near equilib- 
rium (p = 0), the energy increases quadratically as expected. 
The curve exhibits a cusp at 60% shuffle amplitude due to 
a change of magnetism from FM to AFMI. 31 During this 
transition, the magnetic moment — reproduced in the inset of 
Fig. [3] — changes abruptly from 2.2/is/atom. to 1.1/xs/atom 



Energy plots associated with the consecutive mechanism 
are drawn in Fig. EJa). From E Sl (p), one can trivially com- 
pute the energy of shuffle per atom and, thus, can simulate the 
consecutive motion of four layers under the assumption that 
the energy for the shuffling of only one layer is equal to that 
for shuffling all layers, divided by the number of layers. This 
approximation would be correct if there were no interactions 
between layers. For this mechanism, because all atoms are not 
shuffling at the same time, p describes the consecutive shuf- 
fling of all the layers. Formally, from value n — 1 to value n, 
(p — (n — 1)) is, thus, proportional to the displacement of the 
atoms of the n th layer. As a consequence, for p=Q and p=4, 
whatever the mechanism, all atoms are in the same positions. 
The expression for E Gls (p) between n—1 and n is as follows: 

E CfS {p) = EHMp-jn-l))) _ E*( P = A) 

v ' 4 4 

As a consequence, E Gfs (0), E Gfs (1), E Gfs (2), E GfS (3), and 
E Gfs (A) lie on the same line (dashed red in Fig. |4j». Trivially, 
this mechanism reduces the total energy barrier: In Fig. the 
energy barrier E GlSl (p ~ 3.53) -E GlSl (p = 0) is lower than 
the energy barrier E Sl (p ~ 2.30) - E Sl (p = 0). Finally, the 
third curve (E Gl ) corresponds to the explicit calculation of 
the consecutive motion of four layers. 



C. Discussion 

We focus next on the comparison between the model of in- 
dependent layers (CfSi) and the true calculation (Ci). The 
two curves show the same main tendencies with four energy 



4 




Atoms in fixed plane (A - 1 - B) 



Atom in shuffling plane (A - B) 




-A A 



FIG. 5. (Color online) Evolution of local magnetic moment for the 
eight atoms of the supercell as a function of the reaction coordinate 
p. The black, orange, and blue colors approximately depict the bcc, 
hep and mixed bec/hep values for the magnetic moment. The atom 
numbering is related to Fig.|4] Atoms 1, 3, 5, 7 belong to the A layers 
and atoms 2, 4, 6, 8 are successively moving from and A-type layer 
to B-type one. 



FIG. 4. (Color online) (a) Comparison of energies for three differ- 
ent shuffling mechanisms (for a shear of ei=0.5). (Si) is the si- 
multaneous shuffling of all atoms. (CfSi) stands for "consecutive 
from simultaneous", an hypothetical mechanism that corresponds to 
consecutive shuffles within an assumption of independent layers (in- 
teractions neglected), computed with the energies of mechanism Si 
(see text). (Ci) is the true energy computed when layers shuffle in 
a consecutive way (interactions included). Ei and E a represent, re- 
spectively, the interface energy and the activation energy for the first 
shuffle (see text). (C3) is the energy for the consecutive shuffling of 
the atoms for a shear of 63=0.875. (b) Layer stacking during con- 
secutive shuffling: in blue, layers in which the first atomic shell is 
intermediate between the bcc and hep environments. Arrows indi- 
cate the sign of the magnetic moment in each layer, as described in 
Fig.0 



barriers corresponding to the shuffling of the four atoms and 
three metastable structures corresponding to the shuffling of 
only one, two or three layers (p = 1-3). However, the energy 
of the three metastable configurations (for p — 1-3) differs 
by only Ei = E Cl (1) - E cfSl (1) = 60 meV and is, thus, a 
constant independent of the number of shuffled layers. This 
suggests that Ei comes from the two interfaces (in blue in 
Fig. HI between shuffled and unshuffled layers. The number 
of such interfaces is constant and independent of the number 
of layers. 

To assess the electronic origin of these interface energies, 
Fig.[5]represents the evolution of the d-magnetic moment for 
the eight layers in the supercell as a function of p. The transi- 
tion goes along with a progressive change of magnetism in 
which each successive layer shuffle, from bcc- to hep-like 
configuration, involves local jumps from FM to AFMI. In- 
deed, after plane number 2 is translated, there is a local AFMI 
order around this layer. When plane number 4 is translated, 



AFMI order expands over 5 layers, and so forth. Also, the lo- 
cal moments strongly depend on local coordination, in agree- 
ment with the generalized Stoner mechanism for the appear- 
ance of magnetism as a function of bandwidth. 5 For example, 
when the fourth atom goes from a bcc-like local configura- 
tion to a more compact hep-like local configuration, its local 
moment changes from 2.2 pis to -1.1 p,s with no substan- 
tial difference with the simultaneous mechanism. At the same 
time, the moment of the neighboring atom 5 changes from 2.2 
Hb to 1.7 pis only, because atom 5 has now a local environ- 
ment intermediate between the hep and the bcc. Once atom 6 
has moved, atom 5 ends up (for p=3) with a full hep-like local 
environment and a moment of 1.1 ps- This demonstrates: (i) 
the absolute correlation between magnetism and total energy 
in the transition and (ii) the link between interfacial energy 
and layers with intermediate local moment and coordination. 
The interfacial energy only depends on the number of inter- 
faces, which is constant in the present paper. Remarkably, 
energy barriers are not shifted by this interfacial energy. 

For ei = 0.5, one notices that the energy progressively in- 
creases for all the metastable states up to the final energy gain 
of -0.5eV. An additional calculation is performed for a shear 
of £3 = 0.875 (green lines in Fig.|2]and Fig . |4ji . The resulting 
energy exhibits the same characteristics as for t\ = 0.5, and 
notably an interface energy of 100 meV but with an overall 
gain after shuffle of +0.45 eV. For £3 = 0.875 and only around 
p = 3, a new magnetic order — not stable for t\ = 0.5 — 
diminishes the interface energy to 52 meV. It embodies both 
FM and AFM interactions between the layers and highlights 
the energetic proximity of several magnetic orders. Further 
studies, with larger supercells, would be necessary to study 
this magnetic order, but this is not the goal of the present pa- 
per. 
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Two shear thresholds for the bulk transformation can be 
estimated. A simple linear interpolation between the shuffle 
curves Ci and C3 produces the following energy for a shear 
£2 such that t\ < £2 < £3 and the shuffle p: 

£3 - ei 

A lower nucleation-threshold (LT) shear ~ 0.7 is defined 
as the shear above which the gain E(p = 0, €2)—E{p = 4, £2) 
becomes positive. To complete the transformation, the inter- 
face energy has to be overcome for the first layer only, leading 
to shuffling of subsequent layers at no additional cost. For a 
large number of layers, thermodynamical equilibrium, thus, 
drives the transition to completion. Also, an upper-nucleation 
threshold (UT) shear of e% — 0.84 is estimated from the lin- 
ear interpolation. It is defined as the shear above which the 
gain E(p = 0, £2) — E(p = 1,62) for the first metastable 
state becomes positive. Above this shear, even the first shuf- 
fle occurs at no energy cost. Then, all subsequent shuffles 
are strongly favored, possibly leading to propagation and to a 
global transition. 

Consequences on the kinetics of the a-e transition are as 
follows. Shuffle can be viewed as a thermally-activated ki- 
netic process due to vibrations of atoms. 5 ' 7,9 As emphasized 
above, the energy barrier per atom E a , which defines the acti- 
vation energy, is independent of the interfacial energy so that 
the shuffle of a single layer has, thus, no supplementary acti- 
vation cost. In this perspective, the above provides an alter- 
native mechanism for the transition in which the bulk shuffle 
transformation nucleates from one single shuffle event, possi- 
bly initiated by some initial shear £2, such as, e.g., in a shock 
wave. This mechanism considerably reduces the amount of 
energy needed to trigger the transformation, the dimensional- 
ity of the nucleation process being reduced from 3 to 2. 

In order to describe a transition corresponding to a shuffle 



for a domain, a complete thermodynamical description of this 
mechanism would be necessary. This would require relaxing 
atoms after each shuffle, and is beyond the scope of our paper. 

IV. CONCLUSION 

In conclusion, we utilize the fact that shuffle — or opti- 
cal phonon mechanism — and macroscopic compression — 
shear — have in general different timescales, which allows for 
their decoupling. This is true for iron but also for other sys- 
tems such as Zr, Ba, or related transitions, such as Pu^ Sec- 
ond, in iron bcc-hcp transition, we show that the transition- 
path description requires some peculiar magnetic order differ- 
ent from the magnetic order of equilibrium phases. Third, we 
highlight a new mechanism for this transition based on suc- 
cessive shuffle motion of layers. It indicates that each layer 
itself could move by developing "shuffle dislocations" in a 
manner akin to slip motion in plasticity. Our paper suggests 
that detailed investigations of nucleation mechanisms would 
be required to resolve the frontier between the military and 
the thermally activated characters of solid-solid phase transi- 
tions. 
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